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Chimera states in the systems of nonlocally coupled phase oscillators are considered stable in the 
continuous limit of spatially distributed oscillators. However, it is reported that in the numerical 
simulations without taking such limit, chimera states are chaotic transient and hnally collapse into 
the completely synchronous solution. In this Rapid Communication, we numerically study chimera 
states by using the coupling function different from the previous studies and obtain the result that 
chimera states can be stable even without taking the continuous limit, which we call the persistent 
chimera state. 
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The behavior of coupled oscillator systems can describe 
various pattern formations in a wide range of scientific 
fields [11,0. In the systems of nonlocally coupled identi¬ 
cal oscillators, there often appears a strange phenomenon 
called the chimera state, which is characterized by the co¬ 
existence of coherent and incoherent domains, where the 
former domain consists of phase-locked oscillators and 
the latter domain consists of drifting oscillators with spa¬ 
tially changing frequencies 00. This interesting phe¬ 
nomenon was first discovered in the system of nonlocally 
coupled phase oscillators obeying the evolution equation 

^0{x, t) = uj — J dx'G{x — x') V{0{x, t) — 0{x', t)) ( 1 ) 

with 27r-periodic phases 0{x) on a finite interval x G [0,1] 
under the periodic boundary condition, a smooth 2it- 
periodic coupling function T, and the kernel G{y) = 
(/t:/ 2) exp(—/^l^l), where a constant l/n denotes the cou¬ 
pling range j3| . Recently, similar spatiotemporal patterns 
have been found in various systems using, e.g., the logis¬ 
tic maps MM, Rossler systems [HI, and FitzHugh- 
Nagumo oscillators M- 

In the study of the chimera state, the system, Eq. o, 
with the sine coupling [l^l 

r(0) = -sin((/) + (a) (2) 

is particularly important because of its simplicity and 
generality. In fact, this coupling function was used also 
in the first discovery of the chimera state Q . For numer¬ 
ical simulations, we usually discretize Eq. o into such 
a form as Eq. In the simulations of such discretized 
systems, we can confirm that chimera states are surely 
stable in the continuous limit N ^ oc. However, the sta¬ 
bility of chimera states in finitely discretized systems is 
questioned. In fact, it is reported that when N is finite, 
chimera states with the sine coupling are chaotic tran¬ 
sient and finally collapse into the completely synchronous 
solution [T^ . [T^ . . 

Recently, Ashwin and Burylko proposed the weak 
chimera similar to the chimera state, which is de¬ 
fined by the coexistence of frequency-synchronous and 
-asynchronous oscillators in the systems of coupled indis¬ 


tinguishable phase oscillators but is not necessarily spa - 
tially structured as coherent and incoherent domains . 
They studied the weak chimera in some types of net¬ 
works composed of the minimal number of oscillators 
with the Hansel-Mato-Meunier coupling, Eq. 0, and 
demonstrated that the weak chimera can be persistent 
(non transient). In this Rapid Communication, we study 
chimera states in the systems of nonlocally coupled phase 
oscillators with the Hansel-Mato-Meunier coupling by 
numerical simulation, and demonstrate that it is possible 
for persistent chimera states to appear. 

As a model, we consider a ring of N identical nonlocally 
coupled phase oscillators described as 


0j{t) — u; -h 


2R 


( 3 ) 


k=j-R 


with 27r-periodic phases Oj {j = 1,..., A^). This model 
corresponds to a spatially discretized version of Eq. © 
with a constant kernel within a certain range. The natu¬ 
ral frequency uj of the oscillators can be set to zero with¬ 
out loss of generality, and the nonlocal coupling range R 
needs to satisfy 1 < R < {N — l)/2. In this Rapid Com¬ 
munication, we fix ^ 0.35. As the coupling function 
r(0), we choose the Hansel-Mato-Meunier coupling [HI 

r(0) = — sin((/) -P a) Pr sin(2(/)), (4) 

where a is the phase lag parameter of the fundamental 
harmonic component and r is the amplitude ratio of the 
second harmonic component. For r = 0, Eq. © recovers 
the sine coupling, Eq. (j2]). In the systems of globally cou¬ 
pled phase oscillators, it is known that such higher har¬ 
monic components in the coupling function are respon¬ 
sible for a rich variety of synchronous patterns excluded 
by the sine coupling [HhH- Therefore we expect that 
also in the systems of nonlocally coupled phase oscilla¬ 
tors with Eq. 0, we could observe new chimera patterns 
excluded by the sine coupling. 

First, we consider the case of sufficiently large N cor¬ 
responding to the continuous limit. Figure [T] shows the 
results of numerical simulation of Eq. (|3]) with Eq. © 
for several r > 0. For all the simulations of the present 
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FIG. 1. Results of numerical simulation of Eq. © with N = 
2000 and a = 1.46. In each row, the left panel shows the 
snapshot of phase Oj , and the right panel shows the profile of 
the average frequency (Oj) with T = 5000 and trei = 2000 in 
Eq. (|6]). For (a) r — 0.001, (b) r = 0.03, and (c) r — 0.06, 
chimera states are observed, while they are not observed for 
(d) r = 0.12. 


Rapid Communication, we used the fourth-order Runge- 
Kutta method with time interval At = 0.01. In Fig.[Tl we 
fix (a = 1.46, for which chimera states are observed in the 
case of the sine coupling (r = 0). As initial conditions, 
we used 


^i(O) 


6 exp 




( 5 ) 


where Rj G [—1/2,1/2] is a uniform random number, 
which is so close to a chimera state as to assist its emer¬ 
gence [6[. 

In our simulation, chimera states are observed for 
r < 0.073 as shown in Figs. ma)-iii:c). The phase pattern 
(left panels) is clearly separated into coherent and inco¬ 
herent domains, which is characteristic of the chimera 
state. From the right panels of this figure, we can see 
that the average frequency 

1 rtre\-\-T 

{^j)(T) = - ej{t)dt ( 6 ) 

of each oscillator in the coherent domain is almost con¬ 
stant, where T is the measurement time and trei is the 


FIG. 2. Phase diagram of stable solutions to Eq. © with 
Eq. JH at o = 1.46 in the continuous limit {N = 2000). 
Horizontal lines denote the stability regions of each solution. 
We have not clearly determined the transition between the 
chimera state and almost independent oscillation yet. Note 
that the stability region of the wave solution with the wave 
number /c = 1 does not cover r = 0. 


relaxation time, while the frequency in the incoherent 
domain varies continuously. For r > 0.073, chimera 
states gradually disappear as r increases. In addition, 
for r > 0.110, chimera states are not observed, but each 
oscillator evolves almost independently, where the aver¬ 
age frequency seems to converge to a constant value in 
the limit of T ^ oo, though the frequency in Fig. [TJd) 
still exhibits some fluctuations due to a finite T. The 
survey of these behaviors is depicted in Fig. [2l 

From the linear stability analysis, it is found that the 
completely synchronous solution = ^2 = * * * = to 
Eq. (|3]) with Eq. (|4]) is stable for r < (cos(a)/2 0.055 

at a = 1.46). Moreover, chimera states also appear to be 
stable in this parameter region [see Eigs. [T][a) and[T][b)]. 
However, it is reported that when N is finite, chimera 
states at r = 0 are transient and finally collapse into the 
completely synchronous state [T3|. We below confirm 
whether these chimera states, particularly for r > 0, are 
transient or really stable even when N is finite. 

Eigure [3] shows the average lifetime r of the chimera 
state for A = 30, as increasing r from 0 to (cos(a)/2 
0.055. Here we regard the lifetime of the chimera state 
as the time at which the completely synchronous state 
appears, i.e., the global order parameter 


Z{t) = 


1 

N 


N 

k=l 


( 7 ) 


reaches Z{t) = 1. As for the chimera state in the finite 
N cases, it should be noted that it is difficult to judge 
the emergence of the chimera state, because the spatial 
position of the chimera state does not stay still but fluc¬ 
tuates 0, in particular, more violently as N becomes 
smaller. In fact, in the case of A = 30, we could not ob¬ 
serve the characteristic profile of the average frequency 
as in the right panels of Eigs. [Il[a)-[T][c). However, we 
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FIG. 3. Average lifetime r of the chimera state as varying r 
at the parameters A" = 30 and a = 1.46. A point in the fig¬ 
ure is the average over 1000 simulations from different initial 
conditions obeying Eq. ©. 



FIG. 4. Log-log plot of the data (Ar, r) for r* = 0.035 (dia¬ 
mond), 0.039 (circle), 0.045 (cross), 0.050 (square) and 0.055 
(triangle). The data for r* = 0.039 are fitted linearly by the 
least squares method (black line), where we used only the 
data T > 300000 to obtain better linearity for this fitting. 

observed that the coherent domain exists in the phase 
snapshots as in the left panels of Figs. [TJa )-[Hc), which 
convinces us of the emergence of the chimera state. 

For Fig.m it should be noted that there is a possibility 
that chimera states collapse into a stable solution other 
than the completely synchronous state. From the linear 
stability analysis, we can show that the wave solution 
Oi = 0i 27rk{i — l)/N [25| with the wave number 
k = 1 is also stable for r > 0.003 (Fig. [2]). This implies 
that chimera states may collapse into the wave solution. 
However, we never observed such collapse in our simu¬ 
lations from 1000 different initial conditions, Eq. ([5]), at 
each r. 

As r is increased, the average lifetime r increases 
monotonically, and appears to diverge to infinity at a 
certain r = r*. Assuming some values as r*, we ob¬ 
tain Fig. |4]by the log-log plot of the data (Ar, r), where 
Ar = r* — r. From this figure, we can assume the power 
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FIG. 5. (Golor online) Stability region of the chimera state 
(red) on the (coso,r) plane in the continuous limit {N = 
2000). The hatched region corresponds to the persistent 
chimera state, that is, the stability region of the chimera state 
in the case of A = 30. The blue line denotes r = (coso)/2, 
and the completely synchronous solution is stable for r < 
(cos (a)/2. Black circles denote the parameter values of Fig. [T] 
Black triangles denote the parameter values of Fig. [71 where 
multichimera states appear. 

law 

r (X {Ar)~^ (8) 

and we determine r* 0.039 from the best linear 
fitting of the data. Since r* < (cos(a)/2 ^ 0.055, 
this implies that there exists a parameter region where 
the chimera state (with infinite lifetime) and the com¬ 
pletely synchronous state are bistable even in the finite 
A cases. However, we cannot exclude the possibility of 
r* = (cos(a)/2, because it is difficult to obtain the exact 
value of r* due to divergent simulation time. 

Next, we investigate the chimera state of A = 30 for 
r > (cos(a)/2, where the completely synchronous state 
is unstable. The possibility that chimera states appear 
in the region without the stable completely synchronous 
state differs from the case of the sine coupling. In this re¬ 
gion, chimera states cannot collapse into the completely 
synchronous state. Though the wave solution with k = 1 
is stable in this region, we never observed that chimera 
states collapse into the wave solution within our maxi¬ 
mum simulation time t = 2x 10^. Therefore, the collapse 
of the chimera state should not occur if other stable non¬ 
chimera solutions do not exist. Though we searched for 
stable non-chimera solutions other than the wave solution 
by extensive numerical simulations, we could not find any 
such solutions. From the above results, we conclude that, 
in a certain range of r > (cos(a)/2, the chimera state and 
the wave solution are bistable, and the chimera state can 
be persistent (non-transient) even in the finite A cases. 
This result is consistent with r ^ oc for r > r*, as seen 
in Fig. |3l 

Investigating the chimera states in the (cos(a,r) pa- 
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rameter space, we obtained Fig. [5l The red region cor¬ 
responding to the chimera state in the continuous limit 
(A^ = 2000) is spread around r = (cos(a)/2. In the finite 
N cases, the chimera state for small r becomes tran¬ 
sient, while the chimera state for large r remains per¬ 
sistent, as seen at least for r > (cos(a)/2 of the red re¬ 
gion. For cos (a < 0.15, we can see that there exists a 
region r* < r < (cos(a)/2 where the chimera state is 
persistent in the case of = 30. Note that the sta¬ 
bility region of the persistent chimera state (hatched in 
Fig. [5]) extends to the r = 0 line, which implies that the 
chimera state with the sine coupling can be persistent 
(non-transient) even in the finite N cases. Specifically, 
the average lifetime r of the chimera state increases sim¬ 
ilarly to Fig. [Has cos (a is decreased on the r = 0 line, 
and diverge at cos (a* 0.044. However, this fact does 

not contradict the previous study that shows the tran¬ 
sient chimera state [1^ , because the parameter a in that 
study corresponds to the line of black circles in Fig. El 
which has a larger cos a than our hatched region on the 
r = 0 line. 

In summary, we studied chimera states in the systems 
of nonlocally coupled phase oscillators, Eq. m, with 
the Hansel-Mato-Meunier coupling, Eq. by numer¬ 
ical simulations, motivated by the result that chimera 
states with the sine coupling, Eq. ([2]), in finitely dis¬ 
cretized systems are chaotic transient and finally collapse 
into the completely synchronous state [l3| . The existence 
of chimera states was examined in the parameter space 
((a,r) in Eq. (|4]), and the chimera states were observed 
around r = (cos(a)/2 in the continuous limit N ^ oo. 
Eor r < (cos(a)/2, the chimera state and the completely 
synchronous state can be bistable. In this region of the 
finite N cases, the chimera state is transient for r < r*, 
but it is persistent for < r < (cos(a)/2. Moreover, 
even for r > (cosa)/2, it is persistent in the region where 
the chimera state in N ^ oo is stable. At first, we ex¬ 
pected the chimera state to become persistent due to the 
destabilization of the completely synchronous state by 
the effect of r, but have obtained the persistent chimera 
state not only in the unstable region of the completely 
synchronous state as expected but also in its stable re¬ 
gion. As a result, we have discovered that the chimera 
state in the case of the sine coupling can also be persistent 
by using appropriate a in the stability region of the com¬ 
pletely synchronous state. Though we have numerically 
found the persistent chimera state in this Rapid Commu¬ 
nication, its bifurcation-theoretical understanding is still 
an open problem. 

When we investigated the collapse of chimera states at 
a = 1.46, we infrequently observed that a chimera state 
collapses into a weak chimera characterized by the coex¬ 
istence of frequency-synchronous and -asynchronous os¬ 
cillators as shown in Eig.El In [^, the existence 


of weak chimeras for Eq. ([3]) with Eq. (|4]) is confirmed in 
the system with a small number of oscillators {N = 4, 
6, and 10). In our numerical simulation with a larger 
number of oscillators {N = 30), such a weak chimera is 
stable in a small range of r < (cosa)/2, for example, 
0.032 < r < 0.040 at a = 1.46. 

Moreover, as other solutions, we observed multi¬ 
chimera states, which have two or more incoherent do¬ 
mains [H, [23l - [^ , for Eq. (j3j) with Eq. (|4]) in the continu¬ 
ous limit {N = 2000), as shown in Eig.El Other than the 
black triangles in Eig.El we observed multichimera states 
in a large region of the parameter space, though we do 
not describe the region in detail because it is beyond the 
scope of the present Rapid Communication. 
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FIG. 6. Weak chimera for Eq. ([3]) with N = 30, a = 1.46, 
and r = 0.032, which are the parameter values on the line of 
black circles in Fig. El The snapshot of phase Oj (left) and 
the profile of the average frequency (Oj) (right). 



FIG. 7. Multichimera states for Eq. (|3| with N = 2000. The 
left panels show the snapshot of phase Oj , and the right panels 
show the profile of the average frequency (Oj) with T = 5000 
and trei = 2000 in Eq. (|6|). Parameter values are (a) a = 0.95 
and r = 0.28, and (b) a = 1.22 and r = 0.30, which are 
plotted in Fig. El 


[1] Y. Kuramoto, Chemical Oscillation, Waves, and Turbu- [2] A. Pikovsky, M. Rosenblum, and J. Kurths, Synchroniza- 
lence (Springer, Berlin, 1984). 












5 


tion: A Universal Concept in Nonlinear Sciences (Cam¬ 
bridge University Press, Cambridge, 2003). 

[3] Y. Kuramoto and D. Battogtokh, Nonlinear Phenom. 
Complex Syst. 5, 380 (2002). 

[4] D. M. Abrams and S. H. Strogatz, Phys. Rev. Lett. 93, 
174102 (2004). 

[5] S. I. Shima and Y. Kuramoto, Phys. Rev. E 69, 036213 
(2004). 

[6] D. M. Abrams and S. H. Strogatz, Int. J. Bifurcation 
Chaos 16, 21 (2006). 

[7] D. M. Abrams, R. Mirollo, S. H. Strogatz, and D. A. 
Wiley, Phys. Rev. Lett. 101, 084103 (2008). 

[8] E. Ott and T. M. Antonsen, Chaos 18, 037113 (2008). 

[9] E. Ott and T. M. Antonsen, Chaos 19, 023117 (2009). 

[10] E. A. Martens, C. R. Laing, and S. H. Strogatz, Phys. 
Rev. Lett. 104, 044101 (2010). 

[11] O. E. Omel’chenko, M. Wolfrum, and Y. L. Maistrenko, 
Phys. Rev. E 81, 065201(R) (2010). 

[12] M. Wolfrum, O. E. Omel’chenko, S. Yanchuk, and Y. L. 
Maistrenko, Chaos 21 , 013112 (2011). 

[13] M. Wolfrum and O. E. Omel’chenko, Phys. Rev. E 84, 
015201(R) (2011). 

[14] 1. Omelchenko, Y. Maistrenko, P. Hovel, and E. Scholl, 
Phys. Rev. Lett. 106, 234102 (2011). 

[15] 1. Omelchenko, B. Riemenschneider, P. Hovel, 

Y. Maistrenko, and E. Scholl, Phys. Rev. E 85, 

026212 (2012). 

[16] M. R. Tinsley, S. Nkomo, and K. Showalter, Nat. Phys. 
8, 662 (2012). 

[17] A. M. Hagerstrom, T. E. Murphy, R. Roy, P. Hovel, 
1. Omelchenko, and E. Scholl, Nature Phys. 8, 658 
( 2012 ). 


[18] 1. Omelchenko, O. E. Omel’chenko, P. Hovel, and 
E. Scholl, Phys. Rev. Lett. 110, 224101 (2013). 

[19] E. A. Martens, S. Thutupalli, A. Eourriere, and O. Hal- 
latschek, Proc. Natl. Acad. Sci. U.S.A. 110 , 10563 
(2013). 

[20] L. Schmidt, K. Schonleber, K. Krischer, and V. Garcia- 
Morales, Chaos 24, 013102 (2014). 

[21] Y. Zhu, Z. Zheng, and J. Yang, Phys. Rev. E 89, 022914 
(2014). 

[22] S. W. Haugland, L. Schmidt, and K. Krischer, Sci. Rep. 
5, 9883 (2015). 

[23] G. C. Sethia, A. Sen, and E. M. At ay, Phys. Rev. Lett. 
100, 144102 (2008). 

[24] Y. L. Maistrenko, A. Vasylenko, O. Sudakov, 
R. Levchenko, and V. L. Maistrenko, Int. J. Bifurcation 
Chaos 24, 1440014 (2014). 

[25] J. Xie, E. Knobloch, and H.-C. Kao, Phys. Rev. E 90, 
022919 (2014). 

[26] D. P. Rosin, D. Rontani, N. D. Haynes, E. Scholl, and 
D. J. Gauthier, Phys. Rev. E 90, 030902(R) (2014). 

[27] H. Sakaguchi and Y. Kuramoto, Prog. Theor. Phys. 76, 
576 (1986). 

[28] P. Ashwin and O. Burylko, Chaos 25, 013106 (2015). 

[29] D. Hansel, G. Mato, and C. Meunier, Phys. Rev. E 48, 
3470 (1993). 

[30] H. Daido, Prog. Theor. Phys. 88, 1213 (1992). 

[31] K. Okuda, Physica D 63, 424 (1993). 

[32] H. Daido, Physica D 91, 24 (1996). 

[33] H. Kori and Y. Kuramoto, Phys. Rev. E 63, 046214 

( 2001 ). 

[34] P. Ashwin, O. Burylko, and Y. Maistrenko, Physica D 
237, 454466 (2008). 

[35] C. Bick and P. Ashwin, arXiv: 1509.08824 (2015). 



